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with the static hamiltonian and irrelevant operators can be identified 
with the usual 1/M corrections to the static theory. No matching 
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operators with our approach. The values of the matching constants 
come out as a dynamical effect of the renormalization group flow. 
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1 Introduction 

Effective theories for heavy particles are a powerful tool to study heavy fla- 
vor physics. In recent years, many properties of the spectrum and of the 
decays of heavy hadrons have been understood |]l|. The idea underlying the 
effective theory approach is that typical (i.e. after renormalization) momen- 
tum transfers between quarks and gluons in a hadron are of the order of the 
QCD scale Aqcd ~ 300 MeV. As a consequence, relativistic processes for 
a "heavy quark" Q, with a mass M ^ Aqcd, are suppressed. The heavy 
quark is almost at rest and on-shell, and chromomagnetic interactions are 
weak 0. Furthermore, the creation of heavy QQ pairs may be neglected. 

In general, we consider processes where heavy particles are subjected to 
soft interactions. By 'soft' we mean that the energy transfer e and the mo- 
mentum transfer q to the heavy particle H are much less than its mass M: 

I e I, I g I < M. (1) 

Before the scattering H is on-shell and, let us say, at rest. After soft interac- 
tions, H will remain essentially on-shell and will acquire a very small velocity. 
The relevant states for the dynamics of H will be those with 4-momentum 
around the rest momentum 

r = (M,0). (2) 

The effective theories are derived generally by means of an expansion in k/M 
of the lagrangian 0] , where k is the residual momentum defined as 

k = p — r, (3) 

and p is the momentum of the heavy particle. This expansion is of classi- 
cal nature and the effects of quantum fluctuations are taken into account 
imposing a set of matching conditions to the effective theory operators ^ . 

We reformulate the effective theories for heavy particles with a functional 
integral approach 0. We start with a field theory with a cut-off A much 
larger than the heavy particle mass, 

A > M, (4) 

and we integrate away: 



i) the states with a large velocity, for which the spatial momentum | p |~ 
M or greater than that; 

ii) the states which are highly virtual, i.e. the states with an invariant 
mass p^ -C M^ or p^ ^ M^. 

We leave in the effective action Se/f the states with momenta p inside a small 
neighborhood of the on-shell momentum of eq. (^) (see fig.l): 

{po-Mf + p' < A|, (5) 

where Ag is the ultraviolet cut-off of the effective theory, well below the 
heavy particle mass: 

Ae < M. (6) 

In this way we confine ourselves to the low-energy domain. The effective 
action iSeff is defined as: 

exp{tSeff[<l){ky,0<e<Al] }= (7) 

/ n d(j){k) exp{ iS[(f){l)] 0<P < A^] }. 

There exists a different definition of Seff with respect to that one given in 
eq. (^. We may also consider the construction of the effective action for 
heavy particles as a scale transformation. The system is viewed at progres- 
sively larger scales with respect to the Compton wavelength of the heavy 
particle, Xc/'^tt = 1/M. We include therefore in the definition of the effec- 
tive action, in addition to the integration of the high-energy modes, also a 
rescaling of the residual momenta and of the fields. This second formulation 
has a closer relation with the Wilson renormalization group {RG) transfor- 
mation 1^ and consists of the following steps: 

1) We lower the ultraviolet cut-off of the effective action by a factor s, i.e. 
from A^; to Ae/ s. The modes with momenta between A^; and Ae/ s 
are integrated: 

exp{ iSeff[<p{k)- Q<e< [Ae/sY] } = (8) 

/ n d(P{k) exp{ iS[<i){l)] 0<P < A|] }. 

{AB/s)2<fc2<A| 



2) The (residual) moinentum is rescaled according to: 

k' = sk. (9) 

Notice that one does not rescale the energy of the heavy quark, but the 
energy minus the mass. 

3) The heavy particle field is rescaled in such a way that the kinetic 
operator takes a unit coefficient: 

<f>'{k') = C <P{k). (10) 

This transformation can be iterated many times. At each step we lower the 
cut-off A^; by a factor s, i.e. from A^ to Ae/s, from Ae/s to Ae/s"^, etc., 
and we generate a sequence of effective actions Sejf, Sefj, Sefj, ■ ■ ■ ■ 
Expanding the effective action S^ff in a basis of local operators On, 

^eff = Y.^nOn, (H) 

n 

where c„ are given coefficients, the RG transformation has the following 
representation: 

c„ = Jn{Ci , C2 , ..., c^ ,■■■', s) n = 1, 2, 3, ..., k, ... . (12) 

These equations are transformed into differential equations by setting s = 
1 + 6s with 6s <C 1, and expanding up to first order in 6s. 
The fixed point action S* (if any) describes particles which can change their 
momenta with respect to the rest momentum in eq. (^ by a tiny amount 
only, i.e. it describes particles which are essentially static. S* therefore must 
be identified with the static action. A fixed point is a solution of: 

'^n ~ /n(Cx, C2, ..., C^, ...; S) T2 = 1, 2, O, ..., /C, ... . (loj 

Finding a fixed point for our RG transformation can be considered as a 
rigorous proof of the existence of the static limit in quantum field theory. 

Let us make some remarks about the differences between the formulation 
specified by eqs. (p|-|TO|) and the usual one. 

With our formulation we do not deal with lagrangians but directly with 
the functional integral: quantum fiuctuations are taken into account from 



the very beginning. In the usual approach, one performs an expansion in 
1/M of the (relativistic) lagrangian L for a heavy particle to derive the static 
lagrangian Cst- Tree level amplitudes computed with L and with £^4 coincide 
at order (l/M)", implying no need of lowest order matching. Loop diagrams 
are sensitive to high energy physics and turn out to have a different value in 
the full and in the static theory. This implies that quantum fluctuations spoil 
tree level matching, and one is forced to redefine the matching constants in 
order to include their effects. With our approach, there is not any matching 
condition to impose: one has simply to integrate away the off-shell modes 
and look at the modifications of the coefficients of the operators. 

Unlike classical physics or quantum mechanics, the static approximation 
is singular in quantum field theory. The reason is that the heavy particle 
mass M acts as a (physical) ultraviolet cut-off in many processes. Logs of 
M, for example, appear in one loop corrections to heavy currents [Q. In the 
static theory, the heavy particle mass is removed and the resulting divergence 
is regulated by A^;. The cut-off of the effective theory, therefore, does not 
regulate only the divergences related to mass, field and coupling constant 
renormalization, but also those ones induced by the limit M — > cxd. With 
our approach, the singular behaviour of the static approximation is revealed 
by an anomalous scaling as we lower progressively A^;. 

Ordinary RG transformation is a tool to describe sojt interactions of light 
particles in the framework of a field theory. In this case the states with large 
jull energy and with large 3-momenta are integrated away and the effective 
action S'e// describes the dynamics of light particles with momenta inside a 
small neighborhood of the null momentum (see fig. 2) 

n = (0,0). (14) 

We see therefore that our construction of the effective action for heavy par- 
ticles specified by eqs. (|l-0) closely parallels Wilson's RG transformation. 
The difference is that we integrate away the states with momenta far from 
the rest momentum of a heavy particle instead of the null momentum (see 
ref.0 where an intuitive discussion on the above point is presented). 
This paper is organized as follows. 

In Sect. ^ we compute the RG transformation (p|-p!OD in the free case. 
This is the simplest case (exactly solvable) and will serve as an illustration 
of the basic ideas and of the formalism. 



In Sect. ^ we construct the RG transformation for a simple (interacting) 
model, and we discuss the general approximations involved. A perturbative 
evaluation of the RG transformation and of its fixed points is also presented. 

In Sect. ^ we compute the effective hamiltonian as given by eq. (0) - i.e. 
without rescaling - for the same model, and we discuss the relation between 
the original and the effective hamiltonian. A perturbative evaluation of the 
matching constants is presented and the relation with ordinary matching 
theory is treated. 

Sect. ^ contains the conclusions of our analysis. 

2 Free case 

Let us consider the construction of the effective action in a very simple case, 
a free scalar of mass M in Minkowski space. The case of fermions will be 
discussed in a subsequent work. The action is given by: 

'^ = 'C£^^^^P^ [p^-M' + ie] ^{p). (15) 

Since what is small is not the energy of the heavy particle, but the energy 
minus the mass, it is convenient to express the action in terms of the residual 
momentum k. One has: 

where we have defined an 'effective field' d) such that: 



(t){k) = V2M(!>{r + k). (17) 

The difference between the domains of integration of p and k has been ne- 
glected because of the condition stated in eq. (^. In the free case, the 
functional integration is trivial because there is no coupling between high- 
momentum and low-momentum modes, and gives an effective action of the 
same form as the original one, with a smaller cut-off A^;: 

/■Afi d^k 
iSeff = ij^ j—r^ 0HA;) \k^ + e/2M + ie] 0(A;). (18) 

In order to understand the relevance of the operators entering in Seff-, let us 
compute the RG transformation as defined by eqs. (P-p!0|). 
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1) Integrating the modes with momenta between A^; and ^e/ s, the effec- 
tive action becomes: 

iS'^ff = I jy^ -^ <p\k) [ko + ky2M + ze] <P{k). (19) 

plus a constant which does not affect the value of the correlation func- 
tions. 

2) Rescaling the momenta according to eq. @, the effective action reads: 

<ff = ^^-'''^'' r^^'^'^^'/^n^o+^^+^^i^i^v^)- (20) 

3) Rescaling the field in such a way that the kinetic operator takes a 
unit coefficient, 

0'(A;') = s'(^+^)/20(A;7s), (21) 

the effective action looks finally: 

(27r)^^ ^ ^ ^ ° s2M 

We iterate now the transformation many times, say n times; each iteration 
the term with A;^ is multiplied by a factor 1/s and its coefficient reduces 
exponentially with n. There exists therefore a fixed point action S*, which 
is given by: 

/o (2^ 

This is the familiar static action. Note that the mass M of the heavy particle 
disappears in the fixed point action S*. In other words, the actions of heavy 
particles with different masses lie in the same universality class with respect 
to the RG transformation specified by eqs. (p|-p!0|). This is the well-known 
flavor symmetry of the static theory ||^. The term fc^/2M is an irrelevant 
operator with respect to our RG transformation. For high level computations 
this operator may be introduced in the dynamics as an insertion, by means 
of an expansion of the form: 



^S'eff = ^ I 77^<P'\k') W, + -— + ^e] <P'{k'). (22) 



^S* = ^ L J—J^^Kk)[ko + ^e]<P{k). (23) 



giSe// _ giS* 



i+'r(&*'w^^«+ ''<!''"'')]■ p^) 



According to the short-distance expansion (0), the heavy scalar propagator 
is given by: 

^m - ^^ + ^^h^-4)^ + 0(l/M^) (25) 




kQ + ie\ 2MI ko + ie 2M ko + ie 



The round bracket contains a contact term, which is the effect of order 1/M 
of the antiparticle pole, i.e. of the states with a very high virtuality. This 
term can be reabsorbed with a redefinition of the operators containing the 
heavy field 0: 

^-{'"mV (26) 



The second term of the last member in eq. (^5D is the familiar kinetic energy 
operator. We conclude therefore that the irrelevant operators with respect 
to our RG transformation coincide with the operators of order 1/M in the 
expansion of the lagrangian. 

The RG transformation constructed above can be easily extended to de- 
scribe infrared interactions of heavy particles with a fixed velocity f'^ ||^, ^j. 
The only difference is that in the latter case we must leave in the effective 
action the states with momenta in a small neighborhood of the momentum 
Mv instead of the rest momentum of eq. (0) . 

3 A simple model 

We present in this section the construction of the effective hamiltonian in 
an interacting case. We consider a simple model, in euclidean space, whose 
dynamics is determined by the hamiltonian: 

n = d^<l>^d^,<^ + M^<^^<t> + -d^adf,a + -m^a^ + g<l>^<l>a, (27) 

where $ is a heavy scalar field with mass M and a is a light scalar field 
with mass m. This theory is superrenormalizable in four dimensions D = 4, 
because the coupling g has the dimension of a mass, but this does not matter 
for the following considerations. We do not deal with the continuum limit of 



the model (i.e. with the hmit A -^ oo) but rather with its infrared properties. 
This is indeed a model for QED at small energies compared to the electron 
mass. 

We are interested in the soft interactions between these particles; we 
introduce therefore an effective theory where we integrate away the states 
with high-energy for the light particle and the states far from the mass-shell 
state of eq. (||) for the heavy particle. We select a cut-off A^ so that: 



m^ < A| < M^. 



(28) 



Let us compute the RG transformation as specified by eqs. (|^- |10|) . 

1) After making a shift like in eq. ([T6|), the integration of the high energy 
modes is realized according to: 



exp { - H'^ffl <P{k), a{l); < k\ l^ < {Ae/sY] } 



/ 



n d(l){k)da{l) exp { - 

Ae d^l 



As d^k 



(27r) 



D 



4>\k) 



ikA + 5M + 



2M 



<P{k) 



- A 



(27r) 

Ab d^ki d^k2 d^ks 
(27r)^ (27r)^ (27r)^ 



{27rf 6{k, -k2- k,)(j)\k^)<j){k,)a{h)}, (29) 



where A = g/2M is a dimensionless coupling constant in Z^ = 4, and a bare 
mass term 6M has been introduced in the heavy particle lagrangian. 
The integrated modes induce non-local couplings in the effective hamiltonian, 
as well as interactions between an arbitrary number of particles. H'^rr there- 
fore is non local and non polynomial even though the original hamiltonian is 
renormalizable: 



m 



+ 
+ 
+ 



eff - 
^e/s d^l 



Ae/s d^k r k"^ 1 

(2^ '^^^^^ ['^^ + '^^ + 2M " ^^^^^1 '^^^^ 







(2vr) 



D 



a{-l) f + m^-T.L{l) a{l) + 



l'^'^ ^•••^(2^)''^(^i- - ^3) [A - V^^{k,M.h)\<p\k,)<P{k,)a{h) 

^e/s d^ki d^k^ 

(27r)^'"(27r)^ 



{27r)^6{k,... - k,)VZ{ku -, k^)(j)^k,)(j){k2)a{ks)a{k,) 
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+ ... , (30) 

where we have separated the tree level coupling from the loop correction in 
the 00a vertex, i.e. V^^ starts at order A^. 

An exact evaluation of -f^e// is clearly an hopeless task and some ap- 
proximation is needed. Assuming a basis of local operators, this involves a 
truncation of the operator series generated by eq. (pQl). The justification of 
this approximation relies on scale considerations and will be discussed in 2). 

2) Rescaling the full momenta of the light particle and the residual mo- 
menta of the heavy particle according to eq. (P), brings the effective hamil- 
tonian into the form: 



^eff 



^E d^k' 



-(^+l)^t 



(27r 



1 15 



+ 
+ 

+ 

+ 



Af 



Dii 



(27r)^ 
^E d^k'i 

1 , 
^E d^k[ 





■ s' 

/2 



ik', + s6M + 



1 k 



/2 



-sS 



,k\ 



H\ 



s2M • s 

s-(D+2) a{-l'/s) [/'' + s^m^ - s'Si(/7s)] a{l'/s) + 
d^k'. 



f rC ^ 



2n)^'"(2n 



'(27r)^(5(A;;... -4)s"2^ 



D 



kl Kg ^3 ^ 

S 



A - VLC^, -^^ '^ 



5 ) 

s s 



, ru-j . 

' s ■ 



Du 



d^'k 



(27r)^'''(27r) 



D 



{211)'' 5 {k[... 



k'^)s 



-3D 



y y 



Va!f M -^W^)...a^ 



.k[. 

s 



, Ka 



The original limit of integration, A^;, has been restored. 
Let us discuss now the main approximation in the evaluation of H^jf. We 
consider first the free field case (A = 0). As it stems from eq. (0), the mass 
dimensions of the fields 0(/c) and a(/) are given respectively by: 



D + 1 



D + 2 



(32) 



. AjQ 



(31) 



The dimension of a derivative operator is one. Then, the dimension of a 
composite operator O (integrated over space) containing n^ 0-fields , n^ a- 
fields and ng derivatives is given by: 



d. 



o 



D 



D 



D-2 



n^- 



Ua 



rid. 



(33) 



Each RG iteration, the operator O is multiplied by the factor s '^° . There- 
fore, an operator with a positive mass dimension, do > 0, has a weight which 
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reduces exponentially with the number of iterations, and can be neglected 
at leading level. In the case with interactions (A 7^ 0), we assume that their 
effects are small, in the sense that anomalous scaling does not overwhelm the 
canonical one. Also in the latter case, therefore, we can neglect the operators 
with a positive mass dimension. 

The only relevant operators (i.e. with do < 0) are those already present 
in the original hamiltonian. 



0V) (f'^kicj), a^, al'^a, (j)^(j)c 



(34) 



plus 



0^(fc)fci0(A;) 



for 2 = 1,2, 3. 



(35) 



The latter are indeed the operators which enter in the Georgi lagrangian [Q]. 
We will see in section (|3.1|) that they are not induced in -f^e// for the static 
case in a regularization preserving parity. 

At the leading level, therefore, the RG transformation amounts to a renor- 
malization of the original couplings only: 



^eff 



Ae d^k' 



-(^+l)^t 



.k\ 







(27r) 



D 



+ 



+ 



'^ ''' .-(-%f-^l 



(27r)^ 
^E d^k[ d'^k[ 






'dE 



H 



k'^ + s(5M-^H 



dp 



(27r)^'"(27r) 



I^i2rr)^6{k[...-k',)s-'^ 



0/ 
~ A -V- (0,0,0) 



a(-) + 



k' h' 

t(^)0(%a( 
s s 






fvo 



(36) 



3) We rescale now the heavy and the light particle fields in such a way 
that the kinetic operators (p^ {k)ik4^(f){k) and a{—l)l'^a{l) respectively, take unit 
coefficient: 



0'(fc' 
a' {I' 



,-(Z)+l)/2 



-(D+2)/2 



l+l 



'as 



H 



1/2 



dk4 



Hk'/s) 



1- 






0. 

.1/2 



a{l'/s) 



J 



(37) 
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The effective haniiltonian takes finally the form: 



TTlll 



^E d^k' 



A't/;.M 



(27r) 



D 



(j)'\k') ik'^ + 6M' (t)'{k') + 



+ / 7^a'(-n \l'^ + m'\ a\l') + 



(27r) 



D 



(38) 



+ A' 



r^E d^h' rl^h' rl^h' 

, f d k,d k,d k, (2,r)^5(fci-fc^-4)0't(fc;)0'(A:^)a'(A:^), 



(27r)^ (27r)^ (27r) 



_D' 



where the new coupling A' and the new mass parameters 8M' and m^ are 
given by: 



A' 

5M' 



m 



.2-D/2 



1 + Z 



1- 



A 



1 + i 



'dY. 



-\ -1 



/f 



9/C4 






n-i/2 



A 




5M-Sh(0) 
m^ - Sl(0) 



(39) 



Eqs. ( |39D are an approximation of the exact equations (12) for the RG 
flow and, together with (|37D , represent the leading effects of the RG trans- 
formation. To proceed, we must compute explicitly the functions entering in 
eqs. (B^). This is done in the next section with perturbative methods. We 



will search the fixed points of the truncated eqs. (pOf ) in Sect. ^ 



3.1 Perturbative computation 

The construction of the effective hamiltonian presented in Sect. ^ is non- 
perturbative; it is basically a prescription for the degrees of freedom to inte- 
grate away. As such, it does not rely on any Feynman diagram expansion, 
just like Wilson transformation. In deriving it, we assumed only that the 
effects of interactions are not so strong to destroy the distinction between 
relevant and irrelevant operators given by the free theory. This assumption 
is indeed quite reasonable from the physical viewpoint, and it can only be 
checked case by case with an exact computation. However, in order to de- 
velop an intuitive understanding of the properties of our transformation, and 
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to make contact with the usual matching theory, we present in this section 
a perturbative evaluation of the effective hamiltonian. 

Let us sketch the derivation of the Feynman rules for the evaluation of 



H'^rr [|T0|. We decompose the fields and a into their low energy and high 



energy parts: 

0(fc) = (l)Lik) + Mk), ail) = az.(/)+a^(0, (40) 

where the fields with the subscripts 'L' or 'if' contain modes with momenta 
less than or above Ae/s respectively: 

Mk) = ^{k)e{AE/s-k), Mk) = (p{k)9{k-AE/s) 
aiil) = a{l)e{AE/s~l), anil) = a{l) 6(1 - Ae/s). (41) 



We insert then the decomposition (^) in the right hand side of eq. ([29|). In 
the bilinear (free) part of the hamiltonian there are no couplings between 'L' 
and 'ii' fields. On the contrary, the trilinear (interaction) term generates any 
possible couplings between them. We expand the right hand side of eq. (^) 
in powers of A and we compute the functional integral, which is only over 
(pH and an- Each power of 0/, and a^ is associated to an external leg of the 
diagrams, while each bilinear in (pn or a^ is associated to an internal line. 
The external lines of the diagrams have momenta between zero and Ae/s 
while the internal lines have momenta in the range Ae/s — Ae- Loops are 
integrated in a region where all the internal lines have momenta between 
Ae/s and A^. The Feynman rules are: 



— 1 
— : heavy scalar propagator 

k/^ — ik^ /2M 

1 - - - 



(42) 



. light scalar propagator (43) 

-A : vertex. (44) 

Note the asymmetry between the heavy and the light propagator related to 
the shift in the energy for the massive one; the poles for the heavy particle 
and the heavy antiparticle are respectively at: 



k^ = i{MTyk^ + M^). (45) 
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By taking the logarithm on both sides of eq. (pQ]), one sees that only con- 
nected diagrams contribute to H'^ff- Furthermore, one particle reducible di- 
agrams do not contribute to the effective hamiltonian in the low momentum 
region. 

The exact determination oi H'^rr therefore requires to compute and to sum all 
the connected diagrams with an arbitrary number of external legs and with 
generic momenta below Ae/s. Loops induce in the effective hamiltonian in- 
teractions with an arbitrary number of fields and with an arbitrary number 
of derivatives. As explained in the previous section, we neglect higher dimen- 
sion operators since they are irrelevant, and we limit to one loop accuracy. 
Let us go then to the explicit computation of the diagrams. From now on we 
stick to the case D = 4. At one-loop level, the one-particle irreducible dia- 
grams are classified by the number of external heavy and light lines, {uh, ul). 
The self-energy of the heavy scalar [uh = 2, 71^ = 0) is given by: 



T.H{k) 



-tX' 



dH 



D2 (27r)4 U - tP/2M {k - If + w? 



(46) 



where D2 is the region of the /-space, where both propagators have momenta 
between KeI s and A^;: 



{AE/sf < l\ {k - If < A|. 

Since m? = (see later), we can neglect m^ in eq. (0). 
The self-energy at zero external momentum is given by: 



(47) 



Sh(0) 



A^ 



167r2 



:8A, 



1 - - + 0{Ae/M) 



The first derivative of T^H^k) with respect to k^ is given by: 



dki 



k=0 



■2^0 



logs' + 0(As/M) 



(48) 



(49) 



We note that the Taylor expansion of the self-energy around the on-shell 
point A; = is not singular because the lower cut-off Ae/s acts as a regulator 
for the soft divergences. 

First derivatives of T^nik) with respect to spatial momentum components 
ki vanish due to parity. Therefore, linear operators in the spatial derivatives 
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of the form (p^di^) are not induced in the effective hamiltonian by loop ef- 
fects (this was imphcitly assumed in eq. (|36|) , where we have done a Taylor 
expansion of T,H{k) to disentangle the leading terms in H^ff). 
The non-vanishing second derivatives of S// are given by: 



92s 



H 



dkl 



d^T. 



H 



k=0' 



dkf 



k=0 



(50) 



and turn out to be proportional to s/Ae- They are related respectively to 
the irrelevant operators (p'^dfcp and 4>^df(j). 

The self-energy of the light particle Si(/) involves a virtual heavy pair 
and can be neglected at the lowest order in 1/M 



Sl(/) 



0. 



(51) 



The effects of heavy pairs, together with all the other subleading effects [^ 
will be treated in a subsequent work. 
The vertex correction is given by: 



vUkk') 



X' 



dH 



D3 (27r)4 P + m2 
1 



(52) 



1 



k4 + h- i{k + lfl2M ki + U- i{k' + lf/2M'' 



where k and k' denote respectively the momenta of the incoming and out- 
coming heavy particle. At zero external momenta: 



^1(0,0) 



A3 



16772 



21ogs' + 0(A,,/M) 



(53) 



3.2 Fixed points 



Since the first of eqs. (|39|) for the coupling constant A is homogenous, there 
exists a gaussian fixed point: 



A* = 


= 


5M* -- 


= 


2* 


= 0. 



(54) 
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This fixed point is trivial, and we seek a fixed point near this one (i.e. with 
A* -C 1), corresponding to a weakly interacting theory. 
In one-loop approximation, the new coupling constant A' is a function of the 
old coupling constant A only, while the new mass parameters 5M' and m? 
depend on the original mass parameters 5M and m? and on A: 

A' = /(A) 
5M' = g{6M,X) 

m^' = h{m^,\), (55) 

where /, g and h are given functions, which we will compute at one loop 
level later in this section. We note that the equations of Wilson's RG trans- 
formation have the same structure as that in eqs. (^). 
The problem is therefore that of finding a fixed point for the coupling constant 
equation: 

A* = fix*). (56) 

Once this has been found, we tune the original mass parameters in such a 
way that they are unchanged by the transformation; we have simply to solve 
the equations: 

SM* = g{5M\X*) 
m"* = h{m^\\*). (57) 

Unlike the free case, we have 5M* ^ and m?* 7^ 0, in order to cancel the 
contribution to mass renormalization coming from integrated (high energy) 
modes. 

Substituting eqs. (|49|) , (|5l| ) and (^) in the coupling constant equation 
(|39D, we see that the logs^ terms cancel, giving: 

A' = A, (58) 

plus corrections of order K^jM. The coupling constant therefore does not 
change with the RG transformation. This occurs because heavy pairs are 
absent and the model reduces to a free scalar field interacting with fixed 
sources. We have therefore: 







A* = A 

2* n 

m =0 

m* = ^tt^SAe. (59) 
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There exists a fixed point hamiltonian H* with a couphng equal to the orig- 
inal one. Therefore, if we start with an initial coupling Aq <^ 1 at a scale 
~ M, we remain in a region of small coupling, where our perturbative com- 
putation is trustworthy. We see that the formalism developed has a complete 
correspondence with physical intuition. We can consider the result ( ^91) as a 
rigorous proof of the existence of the static limit for our model. 



4 Matching 

We consider now the construction of the effective hamiltonian as specified by 
eq. (^, i.e. without the rescaling of the momenta and of the fields. This is 
indeed the formulation which can be more easily related to the usual approach 
to the effective hamiltonians and to matching theory. We integrate now the 
degrees of freedom with momenta between A and A^, instead of between A^; 
and ^.e/ s as we have done in the previous section. 
At leading level, the effective hamiltonian is given by: 



+ 




(5M-Sh(0) 



0(^) 



(0) 



a{l) + 



(60) 



{2Ti)^5{h - k2 



A -V- (0,0,0) 



<P\k^)<P{k^)a{h 



■(27r)^ 

where a tilde indicates the new range of integration. 

The functions ^//(O), Sl(0), {dtH/dki)o, {dtL/df)^ and V^^iQ, 0, 0) contain 
the leading effects of the integrated modes. We present a one-loop calcula- 
tion of these functions in Sect. 14. 11. Sect. 14. 21 contains a discussion and a 



generalization of the one-loop results. The functions entering in eq. ([60| ) can 
be related to the perturbative matching constants between the full and the 



static theory. This will be done in Sect. O. 



4.1 Perturbative computation 

The Feynman rules are the same as those in Sect. |3.1| with the only difference 
that now loops are integrated in a region where all the internal lines have 
momenta between A and A^;, instead of between A^; and Ae/s. 
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The self-energy of the heavy scalar is given by: 



T^H^k) = —iX' 



dH 



{2TTyh-iP/2M {k-lf + m?' 



(61) 



Let us make some qualitative remarks about the physical meaning of this 
diagram (and of similar ones). Near the upper limit of integration, the loop 
momentum is very large, / ~ A, the shift (^ is irrelevant and the integrand 
is similar to the corresponding one in the full theory. On the contrary, near 
the lower limit of integration, the loop momentum is very small, I ~ Ag. 
We have that /^/2M <C ^4, and the integrand resembles that of the static 
theory. The integrand in eq. ( pTj ) therefore interpolates between the region 
in momentum space in which the heavy particle H is dynamical and the 
region in which H is essentially static. The amplitude therefore includes the 
effects of the fluctuations with momenta both greater and smaller than M. 
The transformation (0) indeed lowers the cut-off in such a way that we cross 
a physical threshold, the heavy particle mass. 
The self-energy at zero external momentum is given by: 



tui^) 



X' 



167r2 



-2M 






4^ + 0((A^/M)2) 



(62) 



The mass renormalization of the heavy scalar 6M is related to the self-energy 
by: 

5M = -S//(0). (63) 

The first derivative of J^h with respect to k^ is given by: 



dE 



H 



dkd 



k=0 



■ A' o 

^ ^2 

167r2 



M2 
log-^-l + 0(A^/M) 



(64) 



In a renormalizable theory (instead of superrenormalizable) , we would have 
an additional term of the form log A^/M^ at the right hand side of eq. (|6^) . 
In our model, a term of this kind appears only in the mass renormalization, 
because field and coupling constant corrections are u.v. finite. The physical 
origin of the two kind of logs will be discussed in sect.^^ 

One of the effects of the integrated modes is that of modifying the nor- 
malization of the kinetic operator given in the original hamiltonian. This 



effect is quantified by the Z-factor relating the on-shell renormahzed field 
0OS to the bare field 0b, defined by 



'^os 



and given by: 



1-i 



Z 



The vertex correction is given by: 

dH 1 



V{k,k') 



X' 



(27r)4 P + m^ 
1 



dkA 



X 



fc=0 



(65) 



(66) 



(67) 



ki + h- i{k + iy/2M ki + h- i{k' + lf/2M 
At zero external momenta: 



mo) 



A3 



167r2 



21og-^ + 0(As/M) 



(68) 



In a renormalizable theory we would have an additional term of the form 
log A^/M^ at the right hand side of eq. (|68D, related to the ultraviolet diver- 
gence of the vertex correction. 



The bare charge is therefore renormahzed as follows: 



A 



1 + 



A2 



M2 

a1 



16,22 log ^ +0(A^/M) 



1 + 



A^ 



167r2 



2 + 0{KeIM) 



(69) 

The cancellation of the log M"^ / K\ terms is related to the conservation of the 
current 0V i^i the effective theory and is not specific of our model. 

In QED, we would have an additional term of the form /5o log A^/M^ at 
the right hand side of eq. (§^). This implies that the (bare) coupling constant 
of the original theory must be evolved from A to M only, because it remains 
unchanged as we enter the effective theory region (i.e. the cut-off becomes 
smaller than M). 
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4.2 Qualitative considerations 

Quite generally, we can identify three regions of momenta which must be 
integrated for the derivation of the effective hamiltonian. 

1) The first one involves momenta / between the original cut-off A down 
to a scale of the order of the heavy particle mass: 

M^ < /2 ^ Al (70) 

In this region, the integrand is similar to that one in the full theory, and the 
mass M of the heavy particle can be neglected. We are basically computing 
the quantum fluctuations for a relativistic massless particle. The contribution 
to mass renormalization coming from this region, for example, is given by: 

5M r dH , A2 

log 77^' (71) 



M J (/2)2 ^M2^ 

because the mass of the heavy particle acts as an infrared regulator. It is 
therefore this region which produces the logA^/M^ term in the mass renor- 
malization constant of our model, and which would give an analogous term in 
Z and V for a renormalizable theory. The physical origin of the logarithm is 
the scale invariance of a massless field theory: there is not any scale entering 
the integrand of eq. (^), so that every order of magnitude range of energies 
gives the same contribution to 6M/M. These kind of logs can be resummed 
with the usual technique of performing many small cut-off lowerings, from A 
to A — 6A, from A — 5A to A — 26A, etc. Region ( [7DD is, let us say, a 'scaling 
region', because we do not cross any physical threshold and we are basically 
in a massless theory. This region is unbounded from above and is the one 
relevant to study the continuum limit of the theory. 

2) The second region is the border between the full and the effective 
theory, and involves momenta of the order of the mass, 

/2 ~ M^. (72) 

This is not a 'scaling' region because we are crossing a physical threshold, 
the heavy particle mass. This region does not produce large logs because it 
does not extend for many orders of magnitude in the energy scale, but gives 
rise to finite terms in the matching constants. Crossing region 2), physics 
goes from a relativistic to a non relativistic one: there is no more enough 
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energy to produce heavy pairs and to give a relativistic momenta to a heavy 
particle. 

3) The third region extends from a scale of the order of the heavy particle 
mass M, down to the ultraviolet cut-off A^; of the effective theory: 

A| < /=^ < M2. (73) 

This region is unbounded from below because A^; may be arbitrarily small 
with respect to M . This is again a 'scaling' region, because the heavy particle 
mass disappears from the dynamics. We are in a different kind of 'massless' 
theory with respect to region 1). The contribution of this region to field and 
vertex renormalization is given by integrals of the form: 



dH , M 



2 



log-T^. (74) 



ll P ^ Ai 

In general, it is this region which produces the logarithms of the form 
log(M^/A|;). The origin of these logs is the scale invariance of a theory de- 
scribing massless particles interacting with fixed (i.e. infinite mass) sources. 
There is not any scale entering in the integrand, and the integral is regulated 
by the heavy particle mass and the effective theory cut-off. These logarithms 
can be resummed with the usual technique, composing small lowerings of the 
cut-off of the effective theory. This region has been extensively studied in 
the previous section. 

Up to now we have considered a model with a heavy particle and a mass- 
less one, but it is not difficult to extend the above considerations to a model 
with, for example, two species of heavy particles and a massless particle. 
Let us assume that the masses of the heavy particles are very different, 
i.e. M' <^ M (this case presents in heavy flavor physics, when both the 
beauty and the charm quarks are treated in the effective theory and their 
ratio Mh/Mc is considered large). In this case, there are two others relevant 
regions to integrate away to derive the effective theory. The first one lies 
between 1) and 2), and consists of momenta between M and M', 

M'^ < /^ < Ml (75) 

In this region the heavier particle behaves as a static particle, while the 
lighter one behaves as a massless particle. This is a 'scaling region'. The 
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contributions to the renormalization constants from this region are therefore 
of the form: 

W ^ °"^' ^^'^ 

because the larger mass M acts as an ultraviolet cut-off to the integral while 
the smaller mass M' acts as an infrared cut-off. A log like that in eq. (|76[) is 
usually called an hybrid log |^. 

The second region is the border between the full theory for the lighter 
particle and the effective one. It consists of momenta / of the order: 

/2 ~ M''. (77) 

We can extend this analysis to a theory with an arbitrary number of 
heavy particles, with masses ...M" -C M' -C M. In every interval between 
two masses, we have a scaling region which gives rise to a log of the mass 
ratio, and a boundary region. 

4.3 Connection with usual matching theory 

Let us discuss now the connection between our formulation of the effective 
hamiltonians and the usual matching theory. 

Consistency requires that the amplitudes of the static theory coincide with 
those of the original (full) theory at the lowest order in 1/M. To this end, 
renormalization constants for the effective theory operators are introduced, 
which are determined by a set of matching conditions. Let us consider a spe- 
cific example, matching of the heavy scalar propagator. The bare propagators 
of the full and the static theory are given, near the mass-shell, respectively 
by: 

— % 

'^^^^^ " A;4[l + z(9S^/9A;4)o]+0(l/M) ^^^^ 



t^s{k) 



-I 



k^{\^%{dY.sldk^)^\-ie- 



where: 




h {2TiY{Pfh-iP/2M ^ ' 



-2tX^ 



^E dH 1 
(27r)4 (/')'■ 
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To avoid cumbersome expressions, mass terms have been omitted in eqs. 
The matching constant Zm is defined requiring that the static propagator 
coincides with the full-one at the lowest order in 1/M: 



i^p{k) = ZmiAsik). 



(80) 



Using eqs. (|78D and eq. (|80D, we derive the following expression for the match- 
ing constant: 

" ^9S^\ fdJ:s\ 



dk/L 



dkA 



(81) 

\'-^"'4/o. 

Let us consider now the renormalization of the heavy particle field induced 
by our transformation. After the functional integration has been done, the 
bare propagator of the heavy scalar is given, near the mass-shell, by: 



iA{k) 



-iZ 



ki-ie + 0{l/M)' 



(82) 



where Z has been defined in eq. (|65|). 

The renormalization constant Z takes into account the modification of the 
field normalization as we scale down the cut-off from A to A^;. This factor 
therefore includes the effects of the fluctuations with momenta between A 
and A^; in the normalization of the fleld, and must be identifled with the 
matching constant Zm'- 

Z = Z^. (83) 



Using eqs. (0), (11), ([79|) and (0), we see that the two matching constants 
do not coincide completely, and the difference is given by: 



Z 




(84) 



iP/2M 



Expanding in powers of 1/M the heavy scalar propagator in eq. (0), we see 
that Zm — Z contains terms of the form {Ae/M)^ <^ 1 with n > 1. Zm and 
Z therefore coincide at the lowest order in 1/M. We believe that this is true 
at every order in perturbation theory. The same result holds for all the other 
Green functions, i.e. differences in the matching constants arise only at order 
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1/M. The static theory is constructed dropping 1/M operators as well as 
corrections of order 1/M to the matching constants of the static operators. 
We conclude therefore that the transformation (|^ reproduces correctly the 
matching constants of the static theory. 

5 Conclusions 

We have formulated the effective theory for heavy particles with renormal- 
ization group techniques, by integrating away the states with high velocity 
and with high virtuality. The fixed point hamiltonian of our transformation 
coincides with the static hamiltonian, and irrelevant operators are identified 
with 1/M corrections. Matching conditions between the full and the static 
theory have not to be imposed by hand as in the usual approach, but come 
out quite naturally as dynamical effects of the RG flow. 

Our interest was to relate, from first principles, quantum field theory, a 
system with an infinite number of degrees of freedom, with quantum me- 
chanics. 

The formulation of the effective theories for heavy particles presented 
is non-perturbative, just like Wilson's transformation. The prescription for 
the degrees of freedom to integrate away and the equations for the RG flow 
indeed do not rely on any perturbative expansion. Our transformation can 
accomodate, for example, instanton effects [|l3l , or can be computed exactly 



with numerical methods. Whenever we apply perturbation theory to it, we 
are able to reproduce all the know results for the matching constants. 

While this approach does not provide us with new computational tools, we 
think that it clarifies the physical content of the effective theories for heavy 
particles, and constitutes a solid framework for (future) non-perturbative 
computations. 
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Figure Captions 

Fig.l: Momentum space of the heavy particle. The states with mo- 
menta inside the white region are integrated away, while the states 
with momenta inside the black sphere are not integrated and consti- 
tute the dynamical variables of the effective action for heavy particles. 
The black sphere is centered around the momentum of a particle at 
rest, r = (M, 0), and has a radius A^ ^ M. 

Fig. 2: Momentum space of the light particle. The states with momenta 
inside the white region are integrated away, while those one inside the 
black sphere are not integrated and constitute the dynamical variables 
of the effective action for light particles. The black sphere is centered 
around the null momentum and has a radius A^; <^ A. 

Fig. 3: One-loop self-energy diagram of the heavy scalar. The contin- 
uous line represents the heavy particle, while the dotted line the light 
scalar. 

Fig. 4: One-loop correction of the heavy-heavy-light vertex. 
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